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ABSTRACT 



The existence of primordial magnetic fields can induce matter perturbations 
with additional power at small scales as compared to the usual ACDM model. 
We study its implication within the context of two-point shear correlation func- 
tion from gravitational lensing. We show that primordial magnetic field can leave 
its imprints on the shear correlation function at angular scales < a few arcmin. 
The results are compared with CFHTLS data, which yields some of the strongest 
known constraints on the parameters (strength and spectral index) of the primor- 
dial magnetic field. We also discuss the possibility of detecting sub-nano Gauss 
fields using future missions such as SNAP. 

Subject headings: Cosmology: primordial magnetic field, weaklensing, structure 
formation 
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1. Introduction 



In recent years, weak gravitational lensing has proved to be one of best probes of the 
matter power spectrum of the universe. In particular, this method can reliably estimate 



the matter power spectrum at small scales which are not directly ac cessible to ot 



her 



methods e.g. galaxy su r veys (for d e tails and further references see e . g. 



Hoekstra fc Jainl (120081 ) 



Refregier 



(120031 ): Bartelmann fc Schneiderl ( I200ll )). 



Munshietal 



( 120081 ): 



Magnetic fields play an important role in the many areas of astrophysics, and are 
ubiquitously seen in the universe. They have been observed in the galaxies and clusters 



of gal axies with the coherence lengths up to ~ 10-100 kpc (for a review see e.g. 



Widrow 



2002) 



T here i s also evidence of coherent magnetic fields up to super-cluster scales 



(IKim et al. 



19891 ) . Still little is known about the origin of cosmic magnetic fields, and 



their role in the evolutionary history of the universe. These fields could have origina ted 



from d ynamo amplification of very tiny seed 



(119791 ): IZeldovich. Ruzmaikin fc Sokolovj ( 119831 ) 



m agnetic fields ~ 10 20 G (e.g 



Parker 



Ruzmaikin. Sokolov fc Shukurov! ( 



19881 )). 



It has been shown that dynamo m echanism can amplify fields to si 



collapsing objects at high redshifts (|Rvu et al 



2009 



de Souza fc Opher 



2010 



Federrath et al. 



2008 



gnificant values in 



2011a 



Sch 



eicher et al. 



2010 



Schober et al 



Arshakian et al. 



201l|). It is also 



possible that much larger primordial 



inflationary phase ([Turner &: Widrow 



magn e tic fie 



1988 



Ratra 



d (~ 10 9 G) were generated during the 



1992] ) and the large scale magnetic field 



observed today are the relics of these fields. In the latter case, of interest to us in this 
paper, magnetic field starts with a large value in the intergalactic medium, while in the 
former case large magnetic fields are confined to bound objects. 

While the presence of primordial magnetic fields have the potential to explain the 
observed magnetic fields coherent at a range of scales in the present universe, such fields also 
leave detectable signatures in important observables at cosmological scales in the universe. 
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The impact of large-scale primordial magnetic field s on CMBR temperature and 



polar i zation anisotropies has been stud i ed in detail ( e.g. 



2002); 



(2005) 



Seshadri fc Subramanianl (12 001); 



Tashiro fc Sugiyama ((2006); 



Sethi & Subramanian (2009) 



(|2Q05|); 



Giovannini fc Kunzd (120081 ) 



Mack et al 



Subramanian fc Barrow ( 



(2002); 



Sethi fc Subramanian (2005) 



Sethi, Hai man. 



andey 



Lewis 



(2010) 



(2004); 



1998B, 



Gopa 



Sethi et al. 



& Sethi 



(2008); 



Kahniashvili &: Ratra 



Yamazakil ( 120081 ); ISeshadri fc Subramanianl (120091 ) ) . 



More recently, lower bounds ~ 10 G on the strength of magnet i c field s 



basec. 



(|2Q10|); 



on observations o f high -e nergy 7-ray 



Tavecchio et al 



bQlpJ); 



p hoton s (e.g. 



Taylor et al. 



km). 



Dolael (12010l ): 



lave been obtained 



Neronov fc Vovk 



Wassermanl (119781 ) showed that primordial magnetic fields can induce density 
perturbations in the post-recombination universe. Further work along these lines have 
investigated the impact of this effect for the formation of first structures, reionization of 



the universe, and the signal from redshifted HI line from the epo c 



Kim et al 



1996 



Gopal fc Sethi 



2003 



Sethi fc Subramanian 



2005 



1 of reionization (e.g . 



Tashiro &; Sugiyama 



2006 



Schliecher. Banerjee. Klessen 



2009 



Sethi fc Subramanian 



2009j). The matter power 



spectrum induced by primordial magnetic fields can dominate the matter power spectrum 
of the standard ACDM model at small scales. Weak gravitational lensing can directly 
probe this difference and therefore reveal the presence of primordial fields or put additional 
constraint on their strength. 

In this paper we attempt to constrain primordial magnetic fields within the framework 
of the two-point shear correlation function induced by gravitational lensing, including the 



contribution of matter perturbations induced 



results with the CFHTLS Wide data (IFu et al. 



)y the se magnetic fields. We compare our 



20081 ). 



Throughout the paper, we used flat (k=0) ACDM universe with Q m = 0.24, O& 
0.044, h = 0.73 and a$ = 0.77. 
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2. Matter Power Spectrum 



Tangled magnetic fields can be characterized by a power-law power spectrum: 
M(k) = Ak n . In the pre-recombination era, the magnetic fi elds are dissipated at scales belo w 



a sca l e corresponding to fcn 



1998 



Subramanian & Barrow 



200 x (10- 9 G/Beff) (e.g. 



Jedamzik. Katalinic. fc Olinto 



1998AI ). Here B e g is the RMS at this cut-off scale for a given 



value of the spectral index, n. Another possible normalization, commonly used in the 
literature, is the value of RMS at k = IMpc -1 , Bq. These two normalizations are related 
as: B cS = -Bo^mi 3 ^ 2 . It is possible to present results using either of the pair {B cSl n} or 
{B ,n}. 

Tangled Magnetic fields induce matter perturbations in the post-recombination era 
which grow by gravitational collapse. The matter power spectrum of these perturbations 
i s given by: P(k) oc k 2n+7 , for n < —1.5, the range of spectral indices we consider here 



( jWasserman 



1978 



Kim et al. 



1996 



Gopal fc Sethi 



20031 ). 



The Magnetic field induced matter power spectrum is cut-off at the magnetic field 



Jeans' wave number: kj ~ 15(10 9 G/B eff ) (e.g. 



Kim et al. 



1996 



Kahniashvili et al. 



2010h . 



The dissipation of tangled magnetic fi eld in the post-recombination era also results in an 



increase in the thermal Jeans' length ( jSethi fc Subramanian 



2005 



Sethi et al. 



2008|). For 

most of the range of magnetic field strengths considered here, the scale corres ponding to k j 



genera lly exceed or are comparable to the thermal Jeans length (Figure 4 of 
(I2008J)). 



Sethi et al. 



For our computation, we need to know the time evolution of the matter power spectrum 
induced by tangled magnetic fields. It can be shown that the domin ant growing mode in 



this case has the same time dependence as the ACDM model (see e.g. 
and references therein) 



Gopal fe Sethi! fl2003f ) 
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3. Weak Lensing & Cosmic Shear 



The cosmic shear power spectrum Pk{t) or the lensing convergence power spectrum, 



P K , is the measure o 



the projection of matter powe r spectrum P$ and is given by the 



following expression (IBartelmann fc Schneider! (120011 )). 



9 



Ho 

c 



\! nn 



dx 



o a 2 (x) Xfxix) 



-;x 



Xhm , , ( , Jk (x -x) 

Jk{x) 



X 



where x is the comoving distance along the light ray and xum is the limiting comoving 
distance of the survey; fxix) is the comoving angular diameter distance; for spatially flat 
(K=0) universe fxix) is numerically equal to the x an d the expression for x m the flat 
universe is as given below, 



= TT !\^m{l + zf + n K )- l / 2 dz (2) 
-"0 Jo 

n(z) is the redshift distribution of the sources and I is the modulus of a two dimensional 
wave vector perpendicular to the line of sight. P$ is the matter power spectrum. In this 
paper, we use tangled magnetic power spectrum as P$ to compute the shear power spectrum 
for the magnetic cases. 

The cosmological shear field induced by density perturbations is a curl-free quantity 
and is donated as an E-type field. One can decompose the observed shear signal into E 
(non-rotational) and B (rotational) components. Detection of non-zero B-modes indicates 
a non-gravitational contribution to the shear field, which might be caused by systematic 
contamination to the lensing signal]*] 

x The presence of primordial magnetic fields will also generate the B-modes of the shear 
power spectrum. Both the vector and tensor modes generated by magnetic fields could 
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These decomposed shear correlation functions can be expressed as: 



where £' is given by 

r°° m ( f r \ 2 \ 

(4) 



f°° di9 ( 

z\o) = z-{o) + 1 (4-12 



and £_ are two-point shear correlation functions which are related to the matter power 
spectrum according to the following relation, 



i±{0) = ^- [°° d£ IP K (£)J 0A (£6) 
^ Jo 



(5) 



6 is the angular separation between the galaxy pairs, and Jo,4 are Bessel functions of the 
first kind. 



4. Shear power spectrum from tangled magnetic field power spectrum 



We use the tangled magnetic field matter power spectrum Pg to compute the shear 
power spectrum P K {t) which in turn is used to calculate £ + , £e and £r using E qs (j3J). 



(jlj) & (JHJ). We have used the same source redshift distribution as in 

lax ^/Ct _|_ 



PuetaL 



(|2008|): 



= A — r ; A 



z b + c 



(6) 



where z r , 



Fu et al. 



6. Values of the parameters a, b, c & A we have taken from the same paper 



(120081 ). Values of these parameters as quoted in the paper are as, a = 0.612 ± 



0.043 ; b = 8.125 ± 0.871 ; c = 0.620 ± 0.065 & A = 1.555. To evaluate the integral © we 



sources these modes. Vector modes are likely to play a more dominant role at angular scales 
of interest to us in the paper. We hope to explore this possibility in a future work. 
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changed the variable from x to z using (jSJ). 



X 



a2 n(z ) 



(7) 



where = £/x(z). again P5 (k,z) can be written as, 



Ps(k,z) = P s (k) x D\z) 



where D(z) is g rowth factor, w hich as noted above is the same as for the flat ACDM mode 



and is given by 



Peeblesl (119931 ): 
5fi 



D z 



[^(1 + zf + n^ 2 



l + z 



We took zii m = 2.5 for our calculations as in 



(fUi + 2) 3 + ao 3/2 



dz 



(9) 



PuetaL 



fl2008f ) 



For comparison, we also compute all the relevant quantities for the linear and non-linear 
ACDM models. For ACDM linear po wer spectrum we used P(k, z) = AkT 2 (k)D 2 (z), where 



the transfer function T(k) is given by 



followed prescription given in IPeacock fc Doddd (119961 ) 



Bond k Efstathioul dl984h . For nonlinear ACDM we 



5. Results 



In Figure 1 we show the tangled magnetic field matter power spectra for a range of 
spectral indices n and magnetic field strengths, B at z = 0. The matter power spectra 
are plotted for k < kj; a sharp cut-off below this scale is assumed for our computation. 
For comparison, we have also displayed the linear and non-linear ACDM matter power 



spectra (the non- 



mear 



Dower spectrum is obtained following the method introduced by 



Peacock fc Doddd (119961 )). The figure shows that the magnetic field induced matter power 



spectra can dominate over the ACDM case at small scales. Possible implications of this 




Fig. 1. — The Matter power spectrum is displayed for the magnetic and non magnetic cases. 
Magnetic field-induced matter power spectra are plotted for k < kj in each case. 



9 



excess have already been studied for 
signal from the epoch of reionization 



early formation of structures, reionization, and the H 



S chli echer. Banerjee. Klessen 



2009 



(Sethi & Subramanian 


2005 


Sethi & Subramanian 


2009; 



Tashiro fc Sugivama 



2006 



Sethi. Haiman. Pandey 



2010) . Here we explore the observational signatures of this excess in the weak lensing data. 



In Figure 2 we show the shear power spectra for magnetic and non-magnetic cases. 
The green and red curves present the shear power spectrum for ACDM linear and nonlinear 
matter power spectra, respectively. The blue curve shows the shear power spectrum for 
the tangled magnetic field power spectrum (B cS = 3.0 nG and n = -2.9). In this figure we 
can see the impact of additional power in the tangled magnetic field-induced matter power 
spectrum as an enhancement in the shear power spectrum on angular scales ~ V . 

The peak of the matter power spectra of both the ACDM model and the magnetic-field 
induced matter power spectra are also seen in the shear power spectra. The ratio of angular 
scales at the peak of the two cases correspond to the ratio of these peaks of the matter 
power spectra: k eq /kj. In the ACDM model the power at small scales falls as k~ 3 , while 
kj imposes a sharp cut-off in the magnetic case. In both the cases, there is power at 
angular scales smaller than the peak of the matter power spectra. But the sharp cut-off in 
the matter power spectrum at k > kj results in a steeper drop in shear power spectra as 
compared to the ACDM case. This cut-off ensures that the magnetic field-induced effects 
dominate the shear power spectrum for only a small range of angular scales. 

In Figure 3, the two-point shear correlation functions and are shown for magnetic 
and non -magnetic cases . As noted in the previous section, we use the parameters of the 



paper of 



FuetaL 



(120081 ) for all our computation, which allows us to directly compare our 



results with their data, shown in Figure 3. 



For detailed comparison with the data of 



FuetaL 



(120081 ). we performed a \ 2 including 



the effect of both the ACDM (non-linear model with the best fit parameters as obtained 



10 



by iFu et al.l ( 120081 )) and the magnetic field induced signal. We fitted the sum of these two 
signals ((£e)b + (£e)acdm) against the CFHTLS data to obtain limits on the magnetic 
field strength Bq and the spectral index n. As seen in Figure 3, the magnetic field induced 
signal dominates the data for only a small range of angular scales below a few arc-minutes. 
However, this can put stringent constraints on the magnetic field model. Our best fit 
values are B = 1.5 nG and n = —2.96. In Figure 4, we show the allowed contours of these 



parameters for a range of Ax 2 



Xt 



Xmin- It should be noted th at B n = is an 



fit to the data because we fix the best fit parameters obtained by 



Fu et al. 



(2008). 



acceptable 
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6. Discussion 



Primordial magnetic fields leave their signatures in a host of observables in the universe. 
Their im pacts on CMBR temperature and polarization anisotropies have been extensively 



studied. lYamazakil (l2010[ ) compute the allowed region in the {B ,n} plane by comparing 
the predictions of primordial magnetic field models with existing CMBR observations. 
Other constraints come from early formati on of structures, Faraday rotation of CMBR 



polar ization (e.g. 



fields 



Kahniashvili et al. 



Schleicher k Miniatil (120111 ). 



20101 ) and reionization in the presence of magnetic 



In addition to the upper bounds on the magnetic field strength obtained by these 
observables, recent results sugges ts that there might be a lower b ound of ~ 10~ 15 G on 



the m a gnetic field strength (e.g 



mm; 



Taylor et al. 



Dolad (120101 ): 



Neronov fc Vovkl (I2010l 



Tavecchio et al 



fl201lh ) . This would suggest that the magnetic field lies in the range 



10 15 < B < a few 10 9 G. This range is still too large for a precise determination of the 
magnetic field strength. 

How do our constraints (Figure 4) compare with the existi n g bou nds on primordial 



Yamazakil (120 lOh ) are stronger than 



magnetic fields? CMBR constraints (e.g. Figure 1 of 
our constraints for n < —2.95. For the entire range of spectral indices above this value, 
we obtain stronger upper limits on B$. Our limits are comparable to bounds obtained 
from the formation of early structures, wh ich also arise from excess p ower in the magnetic 



field-induced matter power spectrum (e.g. 



Kahniashvili et al. 



fj2010h ). 



Can primordial magnetic fields be detected in the Weak lensing data? As seen in 
Figure 3, detection of excess power in the measurement of over what is expected 
for the ACDM model, constrained well from other observations, could be interpreted as 
contribution from primordial magnetic fields. 



The present data is noisy at the scales at which magnetic fields begin to make significant 
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Fig. 4. — The figure shows the allowed region in the (B n : n) p 



((£e)b + (£e)acdm) against the CFHTLS data 



Fu et al 



ane, based on the analysis of 



( 120081 ). The shaded area is the \-o 



allowed region. The three curves (from top to bottom) are contours at 5cr, 3cr and la level. 
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contribution, at least partly owing to errors inherent in ground based measurements of 



shear , e.g. correction for point spread function, etc (e.g. Figure 4 of (ISchrabback et al. 
20101 ); a brief look at this figure might suggest that their measurements would already put 
stronger constraints on magnetic field strength than presented here). Future, proposed 
space missions such as SNAP are likely to greatly improve the errors on thes e measurement s. 



A co mparison of Figure 4 of the white paper on weak lensing with SNAP ([Albert et al. 



20051 ) with the Figure 3 of this paper suggests that SNAP would easily be able to probe 



sub-nano Gauss magnetic fields. 

The magnetic field signal could be degenerate with the overall normaliza tion of the 
ACDM model as measured by erg ; WMAP 7-year data give <7g = 0.801 ±0.030 (( 1L arson et al. 



20111 )). WMAP results are in reasonable agreement with the value of as as inferred by 



the weak lensing data. This error is not sufficient to mimic the much larger signal from 



magnetic field strengths considered in this paper (e.g. Figure 4 of ( ISchrabback et al.ll2010l )). 
However, a more careful analysis will be needed to distinguish the error in erg from the 
sub-nano Gauss magnetic fields. 

One uncertainty in our analysis is that the magnetic Jeans' scale, unlike the thermal 
Jeans' scale which is well defined in linear perturbation theory, is obtained within an 



approximation in whi ch the backreact i on of the magnetic fie 



exactly captured (e.g. 



Kim et al. 



1996 



Sethi fc Subramanian 



d on the matter is not 



20051 ). Even though our 



results capture qualitatively the impact of such a scale, there could be more power on 
sub- Jeans' scale which is lost owing to our approximation of the sharp k-cut off. As noted in 
section 2, the cut-off scale is the larger of the magnetic Jeans' length and the thermal Jeans' 
length. Magnetic field dissipation can raise the temperature of t he medium t o ~ 10 4 K, 



thereby raising thermal Jeans' length of the medium (Figure 4 of ISethi et al. 



(120081 ) for a 



comparison between the two scales for different magnetic field strengths). For B > 10 9 G, 
the magnetic Jeans' scale is the larger of the two scales, as the maximum temperature of 
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the medium reached owing to this process doesn't exceed 10 4 K. In the more general case 
also this would be true as photoionization of the medium by other sources, e.g. the sources 
which could have cause reionization of the universe at z ~ 10, also results in comparable 
temperatures. For magnetic field strengths smaller than considered in the paper, the cut-off 
scale is likely to be determined by thermal Jeans' scale, caused by the photoionization of 
the medium by sources other than the magnetic field dissipation. Our approximation allows 
us to identify important length and angular scales for our study (Figure 2 and 3). However, 
further work along these lines could extend our analysis by taking into account the physical 
effects of sub-magnetic Jeans' scales. 

The analysis of Lyman-a forest in the redshift range 2 < z S, 4 is anoth er powerful 



probe of the matter power spectrum of at small scales (e.g. ICroft et al.l 120021 1 . Primordial 
magnetic field can alter this interpretation in many ways: (a) more small scale power owing 
to magnetic field induced matter power spectrum (Figure 1), (b) dissipation of magnetic 



field can change the therma l state of Lyman-a clouds (e.g. 



Sethi fc Subramanian 



Sethi. Haiman. Pandey 



2010 



20051 ). (c) magnetic Jeans' length can reduce the power at the 



smallest probable scale. We hope to undertake this study in a future work. 
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